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K^^^] in chiral perturbation 



Abstract. We re-investigate the radiative charged kaon decay K 

theory, merging the chiral expansion with Low's theorem. We thoroughly analyze the precision of the 
predicted branching ratio relative to the non-radiative decay channel. Structure dependent terms and their 
impact on differential decay distributions are investigated in detail, and the possibility to see effects of the 
chiral anomaly in this decay channel is emphasized. 



PACS. 13.20.Eb Radiative semileptonic decays of K mesons 
Chiral Lagrangians 



11.30.Rd Chiral symmetries - 12.39. Fe 



1 Introduction 

A particularly useful concept for the investigation of ra- 
diative processes is the decomposition of the transition 
amplitudes in an inner bremsstrahlung (IB) and a struc- 
ture dependent (SD) part. The bremsstrahlung part is 
non- vanishing (in fact, divergent) in the soft photon limit, 
and it can be, according to Low's theorem [1], expressed 
entirely in terms of the corresponding non-radiative am- 
plitude and derivatives thereof; in terms of Feynman dia- 
grams, it corresponds to photon radiation off the external 
charged particles. Only the structure dependent part con- 
tains genuinely new information on the photon coupling 
to intermediate hadronic states. 

Low's theorem was applied to radiative Kg^ decays in 
Refs. [2, 3]. As is often the case for processes where brems- 
strahlung is not forbidden by some mechanism, the IB part 
of the amplitude was found to be largely dominant in the 
partial decay widths, and so several of the earlier studies 
(see e.g. Ref. [4]) rather concentrated on precision tests 
of soft photon theorems. Only with the advent of mod- 
ern high statistics kaon decay experiments has it become 
feasible to measure a relatively rare decay channel to the 
required precision such as to even find effects of structure 
dependent contributions. 

Experiments typically concern the branching ratio rel- 
ative to the corresponding non-radiative process — > 
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where a minimal photon energy E"™' as well as a minimal 

cut 



photon-positron opening angle 9^^ 
frame) are specified. 
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(both in the kaon rest 



Table 1. Experimental values of R for the decay K^^^. A single 
value for the angle cut 0^"' denotes a minimal photon-lepton 
angle, while ranges refer to minimal and maximal values. 



Ref. 


rncut 


/icut 


events 


R X 10^ 


[5] 


lOMeV 


26° - 53° 


192 


0.56 ±0.04 


[6,7] 


lOMeV 


26°- 53° 


82 


0.46 ± 0.08 


[6,7] 


lOMeV 


10° 


82 


1.51 ± 0.25 


[8] 


lOMeV 


[none] 


3852 


1.69 ±0.03 ±0.07 


[8] 


lOMeV 


26°- 53° 


1423 


0.48 ± 0.02 ± 0.03 


[8] 


30 MeV 


20° 




0.63 ± 0.02 ± 0.03 



Experimental results for the relative branching ratios 
R for K^r^^ are shown in Table 1, where we display the 
dependence on the angular cut in degrees, although in 
some of the original publications only the cosine is men- 
tioned. Ref. [5] also includes values for R with different 
photon energy cuts (but always for the same angular in- 
terval) . One should note that the results of the analysis in 
Ref. [6] explicitly refer to the decay probability of an inner 
bremsstrahlung radiative process. We have omitted earlier 
experimental data in Table 1, with statistics determined 
by less than 20 candidate events, see Refs. [9, 10]. 

New experimental efforts that ought to supersede most 
of the previous results are under way at NA48/2 [11] and 
KEK-E470 [12], and ought to record more detailed infor- 
mation than just the partial widths, in particular precise 
photon energy distributions in order to extract informa- 
tion on structure dependent terms. They are likely to sig- 
nificantly surpass the statistics obtained in the most re- 
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cent measurements of the corresponding neutral kaon de- 



cay mode 



[13 15]. therefore more 



precise results on structure dependent terms in K^^^ ought 
to be feasible than obtained for those in -f^es^ in the pio- 
neering investigation of Ref. [13]. 

The appropriate theoretical tool to match such ex- 
perimental refinements is chiral perturbation theory 
(ChPT) [16, 17], the effective field theory of the Standard 
Model at low energies. Radiative /C^s decays were calcu- 
lated up to order p"* in the chiral expansion in Ref. [18] (see 
also Ref. [19] for an earlier tree- level calculation), and in 
Ref. [20] this work was extended and combined with the 
useful aspects of the Low expansion for K^^^^. It is the 
aim of the present article to apply the ideas and meth- 
ods developed in Ref. [20] to -f^Js^ decays. In Sect. 2, we 
present the necessary formalism on the K~^j^^ decay ampli- 
tudes, in particular the IB-SD separation. ChPT results 
on structure dependent terms are discussed anal}ftically 
and numerically in Sect. 3, where in particular we present 
the complete 0(j)^^ results for the axial amplitudes. In 
Sect. 4, we derive our prediction for the ratio of branching 
ratios R and give a detailed account of the uncertainties 
in such a prediction. Sect. 5 discusses the possibility to 
extract information on structure dependent terms from 
differential decay distributions, and in particular on the 
axial anomaly in this decay. Finally, we summarize our 
findings in Sect. 6. 

In this work, we disregard the interesting topic of T- 
odd correlations in K"^^^ decays, which was discussed in 
detail in Refs. [21-23]. 



2 Formalism 

In the following, we consider the decay channel 

- ^\v') e+iPe) vM l{q) [Kt^,] (2.1) 

and its charge conjugate mode. We only study the emis- 
sion of a real photon, that is = 0. 

2.1 The decay amplitude 

The transition matrix element for K'^^^ has the form 
T{K+^) = ^eV:,e'^{qr 

{Vf^iy - Af,^) u{pv) 7" (1 - 75) V{pe) 



+ 



^Peq 



'U'{Pu) 1" (1 - 75) {me- i>e- i) 7/1* APe) 



(2.2) 



The relevant diagrams are displayed in Fig. 1. The first 
term of (2.2) corresponds to diagram a), which includes 
bremsstrahlung off the charged kaon, while the second one 
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a) N^" b) 
Fig. 1. Diagrams describing K'^^^ decay 

corresponds to the radiation off the positron, represented 
by diagram b). We have introduced the hadronic tensors 
and A^^, 



I 

with 



V,A, 



had / 



\K+{p)) 



(2.3) 



(2u7^u - d7^d - sji_,s)/3 



yem 



(2.4) 



(2.5) 



whereas F,, is the K^^ matrix element 

F. = {n'{p')\V^'^{0)\K+{p)) . 
The tensors V^i, and A^^ satisfy the Ward identities 

q'^V^. = F, , qi'A^, = , (2.6) 
which imply gauge invariance of the total amplitude (2.2), 



. 



(2.7) 



The K^^^ amplitude can be decomposed into inner 
bremsstrahlung and structure dependent parts. We re- 
quire these two amplitudes to be separately gauge invari- 
ant, and the SD part to be of order q and higher in an 
expansion in powers of the photon momenta. According 
to Low's theorem [1], the IB terms, which comprise the 
part of the amplitude non- vanishing for small photon mo- 
menta (and in particular the infrared divergent pieces), 
are given entirely in terms of the K^^ form factors /+, /i 
defined by 

F,{t) = ^ (2p'J+{t) + ip- p%h{t)) , (2.8) 

with t = {p — p'Y- [We use the form factors /+, /i instead 
of the more conventional /+,/_= /i — /+ in order to be 
able to write the IB amplitude in a slightly more compact 
form below.] This splitting of the matrix element implies 
a corresponding splitting of the hadronic tensors V^i, and 
Afj,^. The axial correlator A^j^ consists of structure depen- 
dent parts only, and can be expressed in terms of four 
scalar functions Ai, i = 1 . . A, 



e^.p4Aip'fq'' + A2q''W-) 



(2.9) 



A. 
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where W = p — p' — q- (We use the convention £0123 = 
+1.) Note that in comparison to Rcfs. [18,20], we have 
factored out the kaon pole explicitly in the definition of 

the structure fimction ^43. 

The decomposition of the vector correlator reads 



ylB ^ ySD 
y ,,,, ^ ' ixi, , 



where the SD piece is chosen such that 



(2.10) 



(2.11) 



For a given choice of V^^ , the structure dependent part 
of the decay amplitude T in (2.2) is defined to be 



iSD 



Gf_ 
%/2 



X U{p^) {I - -15) V{pe) 



(2.12) 



whereas the brcmsstrahlung part is T^^ = T — T^^ . 

It remains to explicitly construct the decomposition 
(2.10). This is done explicitly in Appendix A, where we 
derive the form of Vj^ in terms of /-|- and /i as 



Pjt{2p'J+{W^) + W^h{W^)) 
pq 



+ ^ (2KA/+ + W^.A/i) + gMt) 



^fi = fi{t) - fi{W^ 



+,1 • 



(2.13) 



A slightly different representation for the IB amplitude 

was already derived in Rcfs. [2,3]. It differs from the one 
given above by terms of order q. An important feature of 
the form (2.13) is the fact that it contains all singularities 
at = in the sense that the complete residue at = 0, 
which is a non- trivial function of the momenta p, p', and 
g, is contained in it, and no terms like (qW)'^ /pq (which is 
formally of order q) are shifted to the structure dependent 
tensor. 

The structure dependent part of the vector correlator 
can also be expressed in terms of four scalar functions Vi, 
« = 1 ... 4, in a basis of gauge invariant tensors according 
to [24] 



V 



SD 



1 
71 



Vi {p'^qi, - p'qguv) + V2 {W/^qt. - qW g,,^) 

+V3{qWp'^W,-p'qW^W,) 



+Vi{qWp'^p',-p'qW^p',) 



(2.14) 



2.2 Kinematics 



We briefly collect the essential information on the kine- 
matics of this decay. The Lorentz invariant amplitudes Ai, 
Vi arc functions of three scalar variables, which we often 
take to be 



s={p' + qf , t = {p- p'f , u = {p-q) 



(2.15) 



These variables are particularly useful in the discussion of 
the analytic properties of Vi, A^. The physical region in 
ife37 decays can be represented as follows: for fixed W"^, 
the variables s, t, and u vary in the range 

<t< {Mk - M^f , 

s±=M2-l(f + M2-M|)(i-VK2) 

±l^X'/\t,M%,M^)X'/\t,0,W'), 
s + t + u = M% + M^ + W^ , (2.16) 
where A is the usual Kallen function, 

A(x, y, z) = x'^ + y"^ + z'^ - 2{xy + xz + yz) . (2.17) 

Note that = M^g denotes the neutral pion mass in 
this work. Varying the invariant mass squared of the 
lepton pair in the interval 



ml<W^< {Mk - M^f 



(2.18) 



generates the region covered by s, t, u in Kf.^^ decays. 
Instead of s, t, u, we also use 

pq/MK = E; , pp'/Mk = El , = {pe+Puf , (2.19) 

where E*, E* are the photon and the pion energy in the 
kaon rest frame. These variables are often useful when 
discussing partial decay widths. 

For the four-body decay K'^^^, one needs five inde- 
pendent variables to describe the kinematics of the decay 
completely. We choose the two additional scalar products 



PPc/Mk =E; , x = Peq/Ml 



(2.20) 



where E^ is the positron energy in the kaon rest frame. 
The dimensionless variable x is related to the angle 91^ 
between the photon and the positron: 



xMl = E; (eI - ^Ef - ml cose%) . (2.21) 

The smallness of the electron mass leads to a near- 
vanishing of x for coUinear electron and photon momenta, 
and hence to a near-singularity, which is avoided by cut- 
ting on the angle 0*^. 

The total decay rate is given by 

r{K+ 7r°e+z/7) = (2.22) 
2MKi2nY /'^LiPs(p;p',Pe,P.,g) El^l' ' 

^ ' spins 

where (iLiPs(p^2''iPe,Piy, 9) is the Lorentz invariant phase 
space element for the K'^^^^ process.^ When performing 



^ For the decay of a particle of momentum p into n particles 
of momenta p\,. . . ,p„, one has 

dLIPs(p;pi,...,Pn) = i'^ip-Y.Pi) n ■ 
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the traces over the spins, wc work with masslcss spinors, 
such that the form factors A3, V3, and f\ drop out, and 
S spins 1^1^ ^ bilinear form of the invariant amplitudes 
Vi, Ai, i = 1, 2, 4:, and /+. The explicit result is displayed 
in Appendix B. 



3 Structure dependent terms in ChPT 
3.1 Analytical results at order 

In Ref. [18], the chiral expansion was carried out to or- 
der p"^ for both the neutral and the charged decay modes 
(one-loop order). We do not describe that calculation in 
any detail and only quote the result, adjusted to the sep- 
aration between bremsstrahlung and structure dependent 
terms. The axial amplitudes are given in terms of the 
Wess-Zumino-Witten anomaly term [25, 26], the result is 



-AAn 



A, = 



1 



27r2F2 ' 



A, = . (3.1) 



The vector structure functions Vi were shown to be given 
as 



V, = V2I2 , V2 = -^{h +P'(lh) , 



(3.2) 



where the fimctions .f2(W^) are defined in Ref. [18]. 

It was pointed out in Ref. [18] that these functions 
are real throughout the physical region. In fact, cuts only 
start at t = (MK+M^f, W"^ = (Mk+M^)'^, therefore far 
from the phase space boundaries, which is why the Vi at 
this order are smooth and well-behaved, and even constant 
to a high degree of accuracy. They can therefore be ap- 
proximated by a seemingly drastic simplification, namely 
by their values at the special kinematical point s — M^, 
u = -M|-, t = = 0. This corresponds to an expansion 
to leading order in the photon momentum q, plus setting 
t = 0. The rather compact and simple result reads 



Vi 



(1 -a;)-2 n 



F^- 32.2^. {J (53 -25. + 2.^) 

log a; 
(1^ 



9 -^x log a; 



fl27- 93a; + 21.2 



V2 = -j^{U + L^o) 

{i + x){i-xy 



(l-a;)-3 r 166 



1 + x 



2(1 -x) 

2a; log X 



0{q,t) 



1-x 



loe 



+ x(3 + 2a;)f^-9(12-a;)(4-a;)2 



1-x 



1-x 



+ 0{q,t) , 

+ a;(2 + 3a; + a;^)^^ - 27(7 - a;)(4 - xf ''^^ 



1-x 



+ 0{q,t) 



1-x 
(3.3) 



where x = M^/Mj., y = /M]^, and we have made fre- 
quent use of the Gell-Mann-Okubo relation. See (D.6) for 
a definition of the scale-independent low-energy constants 
ig, Lio- We remark that the expressions for Vi and V3 in 
(3.3) are identical to the equivalent approximations in the 
neutral kaon decay mode quoted in Ref. [20]. 



3.2 Analytical results at order 

As the above-described results for the structure func- 
tions Vi, Ai at 0{p'^) are the leading contributions in 
the chiral expansion, and as the chiral expansion involv- 
ing strangeness does not necessarily converge very fast, it 
is mandatory to understand the structure of higher-order 
corrections in order to give realistic estimates of the struc- 
ture dependent terms. 

The analytic structure of the tensors V^v, A^^ was 
investigated in detail in Refs. [20, 23]. The main results of 
those investigations are: 

1. Due to strangeness conservation, cuts in the variables 
t, u, can all start only at {Mk + M^)"^, which is far 
outside the physical region. As in the case of the vector 
structure functions at 0{p^), these cuts are expected 
to hardly affect the momentum dependence of the Vi, 
Ai inside the decay region, where they can therefore 
be approximated by polynomials. 

2. At higher order in the chiral expansion, there are two - 
and three-pion cuts in the s-channel that make the 
structure functions complex. These cuts were discussed 
in great detail in Ref. [23] . Due to the photon in the fi- 
nal state, all the diagrams developing imaginary parts 
have a P-wave characteristic, i.e. the imaginary parts 
rise only slowly above threshold, and a cusp-like struc- 
ture in the real parts is smoothed out (oc (s — AM^)^^'^ 
in the case of the two pion cuts). For our purposes, 
the real parts of the structure functions can therefore 
still be regarded as "smooth" . 

The main purpose of the analysis of higher-order correc- 
tions is therefore to see how large the corrections to the 
averaged (constant) structure functions might be. 



3.2.1 Complete order corrections to the axial amplitudes 

We have calculated the complete 0{p^) corrections to the 
axial structure functions Ai, A2, and A^^. The contribution 
of A3 to the squared matrix element is always suppressed 
by a factor of m1/M]^ « 10~^ and is therefore neglected. 
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The generic structure is as follows: 



The explicit forms for the various loop functions as well as 
the combinations of low-energy constants entering the ex- 
pressions (3.4) can be found in Appendix C. In particular, 
the constant terms Xi^2 in (3.4) contain terms propor- 
tional to the strange quark mass, and are therefore po- 
tentially big. We remark that by rcnormalizing the order 
amplitudes according to Fj^Fk, all dependence 

on the low-energy constants L4 and L5 is absorbed in the 
physical meson decay constants. 



Table 2. Values for the structure dependent terms as given 
by C(p^) and estimated from 0(/) ChPT. The symbols (Vi), 
(Ai) denote the averages of the real parts of the Vi, Ai over 
phase space, in units of A^k- For the Vi, we also show the ap- 
proximation given in (3.3). For the error bands, see discussions 
in main text. 







(3.3) 


0(p«) 




-1.24 ± 0.004 


-1.23 


-1.24 ± 0.4 


{V2} 


-0.19 ± 0.007 


-0.21 


-0.19 ± 0.2 


(Vs) 


-0.02 ± 0.001 


-0.02 


± 0.1 


{V4) 








± 0.1 


(Ai) 


-1.19 




-1.29 ± 0.4 


{A2) 


0.30 




0.33 ± 0.1 


{As) 


-1.19 






(A,) 







± 0.3 



3.2.2 Polynomial corrections to the vector amplitudes 

A complete evaluation of the vector structure functions 
at order requires a full two-loop calculation and is be- 
yond the scope of this artick\ In order to gain a basic 
idea about potential higher-order corrections, we content 
ourselves with an investigation of just polynomial contri- 
butions. Note that the loop contributions with cuts in the 
physical region are known to be tiny at this order [23]. 

The polynomial part for the Vi can be calculated from 
the Lagrangian Cq [27,28]. For the form factors Vi, V2, 
it contains all possible terms linear in s, t, u, W^, M^, 
M^. The numerically potentially largest corrections for 
Vi , V2 turn out to be the terms suppressed by a factor 
M'^/{A-kF)'^ with respect to the leading Lg, L\q contribu- 
tions. The (leading) polynomial contributions for V3, V4, 
which only appear at 0{p^), are constant. 

Wo have furthermore calculated the contributions of 
the form Lj x Lj at order p^ . In analogy to what was 
done above for the axial structure functions, we find that 
a renormalization of the order couplings according to 
F"^ F^Fk takes care of all such terms. This informa- 
tion is used for the numerical evaluation of the structure 
functions below, where such a normalization is chosen in 
order to minimize higher-order corrections. 

3.3 Numerical evaluation of structure dependent terms 

In Table 2, we summarize the numerical evaluation of the 
various structure dependent terms. We use the parameters 
and numerical values specified in Appendix D. In partic- 
ular, we neglect any variation in the order p'^ low-energy 
constants ig, Liq, the effect of which should bo generously 
covered by the uncertainty in the order p^ contributions 
estimated below. 

The averages (Vi), (Ai) in the first column of Table 2, 
referring to the values of the structure dependent terms at 
order p^, are obtained by integrating the structure func- 
tions over phase space and dividing by the phase space 



voliun(^: the quoted uncertainties are the Icr errors of this 
averaging procedure. These are therefore no realistic es- 
timates of the uncertainties of the mean values, but only 
serve as an illustration to what extent the approximation 
of the Vi being constant at this order is justified. It is 
seen that the variation within physical phase space is ab- 
solutely negligible. This is, on the one hand, due to the 
absence of cuts, but also to the fact that Vi and V2 are 
(at the scale of the p mass) numerically dominated by 
counterterm contributions, which are in turn necessarily 
constant at this order. For comparison, we also show the 
values of the Vi at the special kinematical point s = M^, 
u = M|-, t = W"^ = as given in (3.3), which are consis- 
tent with this picture. 

A numerical assessment of the corrections at next-to- 
leading order (p®) is much more difficult, essentially due 
to the large number of unknown low-energy constants. 
For the axial structure functions that we have calculated 
completely, we proceed as follows: by averaging the (real 
parts of) the loop contributions at a scale n = Mp over 
phase space, we calculate the shifted mean values (Ai); 
the counterterm contributions are added as the essential 
uncertainty. This uncertainty is determined by examining 
the scale dependence of the combined counterterm contri- 
butions. If the scale is varied such that the corresponding 
logarithms change by one, we find for the counterterm 
parts of ^1,2 

^i.ct = T i92^/j,2j,^ - M^) +4s + t + u} , 

As there are no loop contributions to A4, the correspond- 
ing coimterterm combination is separately scale indepen- 
dent, and we have to utilize an even simpler order-of- 
magnitude estimate with the result 

16 
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In order to obtain simple error ranges for (Ai 2), we again 
average the momentum dependent terms over phase space. 
It is seen that the largest contributions to the uncertainty 
stem from constant terms oc Af^. In this sense, within the 
accuracy at which these structure functions can presently 
be predicted, we can even neglect the momentum depen- 
dence and approximate the structure functions by con- 
stants. 

The results thus obtained for the Ai are displayed in 
the third column of Table 2. We note that the uncertain- 
ties for A\^2 are of the size of typical chiral SU(3) correc- 
tions of about 30%. 

For the V^, we use the same order-of- magnitude 
arguments concerning higher-order contributions as in 
Ref . [20] . The uncertainty for the dominant structure func- 
tion V\ is estimated to be of the order of 30%; as V2 is 
suppressed at leading order, we scale its uncertainty by 
a factor of 2. The (constant) counterterm contributions 
to Vz^i are estimated by dimensional arguments similar to 
that for discussed above. All the numbers are collected 
in the third column of Table 2. 



Table 3. Coefficients for the K'^^ phase space integral. 





01 


02 


as 


a4 


0.09653 


0.3337 


0.4618 


3.189 


6.278 



with re = ml/Mj., = M^/Mj.. The expression for TZ 

in terms of these reduced phase space integrals is then of 
the same form as for the neutral decay mode [20], 



■R = 



8a J cJlips SM 
^ JdydzA{y,z)f+{ty 



(4.5) 



We will occasionally also refer to a ratio TZ^^, which is 
understood to be calculated according to (4.5), with all 
structure dependent contributions in SM omitted. 



4.1 Phase space integrals 

Assuming^ 



4 The ratio R 

The ratio R defined in (1.1) is a particular useful quantity 
to consider in K^a-y decays, as it is both the quantity that 
is naturally measured in experiment (as opposed to the 
branching ratio r{Ke3j)/ Fi^ii), and can be predicted in a 
clean way theoretically. 

We repeat here the corresponding discussion of R for 
the neutral kaon decay in Ref. [20] for the decay channel 
K't^^. As the layout of the formalism is analogous to the 
neutral channel, we shall be rather brief and only comment 
in more detail on the numerical results. 

For the moment, we neglect radiative corrections and 
isospin breaking and denote R in the absence of virtual 
and real photon corrections by TZ, 

a-'n=[a-'R]^^^. (4.1) 

We will comment on radiative corrections in Sect. 4.4. 
We define the quantity SM by 

r(K^)^ ''"\^r'-^' Mofl^..sm. (4.) 

such that / dLiPsSM is dimensionless and contains no fur- 
ther (electroweak) coupling constants. By factoring out 
/+(0)^, SM only contains the normalized form factor 
f+{t) = f+{t)/f+{0). The non-radiative width FiK^^) is 
conventionally written as 



7r=^ 



(4.6) 



riK+) = J dydzp{y,z) 



where y = 2ppe/M^, z = 2pp'/Ml^, and 

A{y, z) = 4{z + y- 1)(1 - y) + re{4y + 3^-3) 
- 4r„ + re{r^ - re) , 



(4.3) 



one may expand the integral in the denominator according 

to 

I = j dydzA{y,z)U{tf (4.7) 

= ao + aiA+ + ai (A^ -h 2A'|) -|- asA+A" a4A'|^ , 

where the numerical values for the Uj are collected in 
Table 3. We point out that we have used the physical 
(charged) kaon and (neutral) pion masses everywhere, 
such that the kinematics in the phase space integral cor- 
respond to the physical situation. 

Similarly we can expand the phase space integral for 
the radiative decay in terms of Kg^ form factor parameters 
according to 
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^ duPS SM (4.8) 

= bo + hX+ + hXl + baX'l + hX+X'i + hX'l^ . 

The integral /'*' and the coefficients 6j depend on the ex- 
perimental cuts ii"™', 6'™*. The numerical results for the 
standard cuts £;™* = 30 MeV, 61^!;' 20° are displayed in 
Table 4. Where applicable, the coefficients have been de- 
composed into their bremsstrahlung and their structure 
dependent parts. The relative size of the structure depen- 
dent contributions as predicted by ChPT is very similar 
to the neutral kaon decay channel, they reduce the width 
by about 1%. The uncertainties quoted in Table 4 refer to 
the estimated higher-order contributions in the structure 
dependent terms as discussed in Sect. 3.3. 



(4.4) 



Note that, while wc use M-,, = M^o elsewhere in this text, 
the expansion of the form factor /_|_ is conventionally normal- 
ized to the charged pion mass M^± . 
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Table 4. Coefficients for the K'^^-y phase space integral. 



bo Oi 


62 bf 


64 65 


1.019 3.98 


5.81 11.83 


41.9 84.8 


.SD 


6?° 


bT 


-0.012 ± 0.004 


-0.03 ±0.01 


-0.10 ± 0.03 



Table 5. Coefficients for the A+, A+ dependence of 72. ,72.. 
The error margins for the 72 coefficients are due to uncertainties 
in higher-order contributions to the structure dependent terms. 



72™ (1,0) ■ 10^ 


• 10^ 


cP ■ 10* 


cf ■ 10* 


0.640 


12.0 


-5.4 


16.6 


72(1, 0) • 10^ 


ci • 10^ 


C2 • 10* 


C3 • 10* 


0.633 ± 0.002 


12.5 ±0.4 


-5.4 ±0.3 


16.9 ± 0.4 



4.2 Form factor dependence of 7?. 

We can now study the dependence of TZ {Ti}^) on the form 
factor parameters A+ — A+/A^j_, A" = A"/(A^j_)^, where 
we choose A^j_ = 0.0294 as a central value for the slope 
parameter A+. We expand TZ according to 

n{\+,\'\) = 7^(l,o){l + Cl (A+-i)+c2 (A+-1)' 

+ C3A'; + ...}, (4.9) 

and V}^ accordingly (with expansion coefficients cP). We 
show the numbers for the coefficients in Table 5. 

We find that although there is a significant cancellation 
between the A+, A" dependence of numerator and denom- 
inator, the cancellation is not quite as complete as for the 
-^e3(7) channel [20]. Going from a point-like form factor 
(A+ = 0) to the physical one A+ = A^ induces a change 
of about 1.3% in 7^ or TZ^^ . 

On the other hand, TZ and 7^^^ are remarkably stable 
within the range of uncertainties of the measured slope pa- 
rameters. Even if the latest experimental results that de- 
termine A" [29-32] do not agree perfectly with each other, 
all these results combined show a strong (anti)corrclation 
between A+ and A^ (see e.g. Fig. 8 in Ref. [32]). If we vary 
(A+, A" ) in the parameter space indicated by the different 
IfT ellipses of Refs. [29-32], we find that TZ differs from 
TZ{1,Q) by less than a pcrmille. With the present knowl- 
edge of the form factor f+{t), it is therefore already pos- 
sible to predict TZ to excellent precision, and the biggest 
uncertainty (of the order of 0.4%) stems from unknown 
higher-order corrections in the structure-dependent terms. 

4.3 Dependence on the experimental cuts 

We briefly study the dependence of the parameters in (4.9) 
on the experimental cuts by displaying their values for the 



alternative cuts S™' = lOMeV, 6^]^ = 10° in Table 6. In 
addition, for historical reasons we also show results for 
the angle range 26° < 6**^ < 53° in combination with 
£;™' = lOMeV, compare Table 1. 7^, TZ^^ of course de- 
pend strongly on the cut values. The dependence on the 
form factor parameters, however, remains very mild also 
for these different cuts. In particular, the ratio TZ is even 
more stable with respect to variations of the form factor 
parameters, i.e. the variation of A+, A" in the range of the 
latest experimental results, as described above, leads to a 
change in TZ below the permille level. As the difference be- 
tween the coefficients Ci and is as small for all cuts as 
for the standard ones discussed in the previous subsection, 
we refrain from showing the c-^ in Table 6. 

4.4 Isospin breaking and radiative corrections 

As we can predict the ratio TZ to surprising precision of 
about half a percent, using Low's theorem, the experi- 
mentally available information on the Kg^ form factor /_|_ , 
and ChPT for the structure dependent terms, we have 
to comment on isospin breaking corrections, generated by 
real and virtual photons and the light quark mass differ- 
ence niu — md ^ 0, which may clearly have an impact at 
the percent level. 

As radiative corrections comprise the radiation of (ad- 
ditional) soft real photons, we have to specify the precise 
meaning of the ratio R as defined in (1.1), in the presence 
of virtual and real photons. We define the denominator 
to be the fully inclusive width — > i:^e'^V(.(n'y), where 
(717) denotes any number of photons of arbitrary energy. 
In analogy, the numerator in (1.1) refers to the measure- 
ment of radiative iiT^ decays with at least one photon 
fulfilling the cut requirements E* > E^™', 6**,, > plus 
arbitrary additional photons. 

A complete calculation of radiative corrections in R is 
beyond the scope of this article. Below, we argue that the 
most sizeable effects can easily be taken care of, and that 
the unknown corrections can reasonably be expected to 
be small. 

1. As was already argued in Ref. [20], the large short 
distance enhancement factor Sew oc log Mz/Mp that 
renormalized the Fermi coupling constant Gp [33-35] 
applies identically in numerator and denominator of 
R, so it cancels in the ratio. 

2. In Ref. [34], isospin breaking corrections in were 
calculated in chiral perturbation theory up-to-and- 
including O {p^ , {tUu — md)). Parts of these can 
be collected in a modified value for /+(0), which still 
cancels in R. In addition, the parameters describ- 
ing the slope parameter expansion of the K^s phase 
space integral are shifted, the modified coefficients of 
the expanded phase space integral are given in Ta- 
ble 7. As in Ref. [34], the contribution of (single) real 
photon radiation was restricted to pion/electron mo- 
menta in agreement with non-radiative kinematics, we 
re-modify the coefficient ao to be fully inclusive again, 
increasing it by 0.51%. 



8 



B. Kubis, E. H. Miiller, J. Gasser, and M. Schmid: Aspects of radiative K'^^ decays 



Table 6. 72. ,72. for different values of the experimental cuts on iJ!^"*, Se"', as well as coefficients for the A+, A" dependence 
of 72.. The error margins for the 72 coefficients axe due to uncertainties in higher-order contributions to the structure dependent 
terms. 



^cut 


zjcut 


72™ • 10^ 


7^• 10^ 


ci • 10=^ 


C2 • 10-* 


C3 • 10^ 


30MeV 


20° 


0.640 


0.633 ± 0.002 


12.5 ± 0.4 


-5.4 ±0.3 


16.9 ±0.4 


SOMeV 


10° 


0.925 


0.918 ±0.002 


11.1 ±0.3 


-4.7 ±0.2 


15.0 ±0.3 


lOMeV 


20° 


1.211 


1.204 ± 0.002 


7.5 ±0.2 


-3.2 ±0.2 


10.1 ±0.2 


lOMeV 


10° 


1.792 


1.785 ± 0.002 


6.7 ±0.2 


-2.8 ±0.1 


9.0 ±0.1 


lOMeV 


26° - 53° 


0.554 


0.553 ± 0.001 


5.7 ±0.1 


-2.4 ±0.1 


7.5 ±0.1 



Table 7. Coefficients for the K'^^ phase space integral, includ- 
ing corrections of 0{a, rUn —md)- The values for ai,2 are taken 
from Ref. [34]. We are grateful to V. Cirigliano for providing 
us with the numbers for 03,4, which are not included in that 
reference. For ao, see text. 



oo 


ai 


a2 


as 


04 


0.09583 


0.3287 


0.4535 


3.124 


6.136 



3. The part of this estimate that is necessarily incom- 
plete concerns isospin breaking corrections to the nu- 
merator of R. Radiative corrections that are poten- 
tially large are those enhanced by electron mass sin- 
gularities. As we have defined the numerator of R to 
be inclusive with respect to additional photon radia- 
tion, such electron mass singularities should, accord- 
ing to the KLN theorem [36,37], be taken care of by 
evaluating the running coupling constant at the scale 
of the kaon mass, a a(l + \og{M'^ / . In 
Ref. [34] , a part of these large logarithms are absorbed 
into /+(0), which reduces the correction in the numer- 
ator to a ^ a{l + ^ log(Af |./m^)). 

4. We expect the remaining, non-enhanced, radiative cor- 
rections to be small, of the order of qj/tt. As a conser- 
vative estimate, we assign a relative uncertainty due 
to these of the size ±Ziem = ±5a/7r w ±0.01. 

5. Finally, corrections in the numerator due to the light 
quark mass difference that are not local and absorbed 
in /+(0) are expected to be tiny and are neglected. 

Summing everything up, we find that the corrections in 
the denominator of R plus the remaining electron mass 
logarithms enhance R by 1%. The estimated uncertainty 
in other radiative corrections is larger than that from 
higher-order chiral corrections in the structure dependent 
terms. Our combined predictions, for the various cut com- 
binations discussed before, are summarized in Table 8. 

We briefly compare the predictions in Table 8 to the 
experimental results shown in Table 1. We find very good 
agreement with Ref. [5] as well as with the standard cuts 
result in Ref. [8], while the values for R reported in Ref. [6, 
7] as well as the one for the angular cut range 26° — 53° 
in Ref. [8] are a bit below the theoretical values, though 
disagreement hardly exceeds la deviation. Clearly, forth- 



Table 8. Combined results for R for various cut combinations. 



^cut 


/icut 


R-IO'^ 


30MeV 


20° 


0.640 ± 0.008 


30MeV 


10° 


0.928 ±0.011 


lOMeV 


20° 


1.217 ±0.014 


lOMeV 


10° 


1.804 ± 0.021 


lOMeV 


26° - 53° 


0.559 ± 0.006 



coming even more precise determinations of i? [11, 12] are 
very welcome. 

In Ref. [8], the ratio R is also given with no cut on 
01^, see Table 1. It is difficult to provide a precise theo- 
retical value of R for this situation, because the electron 
mass singularity renders the numerical integration over 
phase space unstable. While our calculation suggests a 
rather bigger value for R than what is quoted in Table 1, 
we refrain from working out a precise result. We do not, 
therefore, provide a final number in Table 8. 

5 Structure dependent terms 
in differential rates 

5.1 Photon energy distribution 

It is obvious from the results of the previous section that 
a precise measurement of the ratio R is insufficient to de- 
termine structure dependent terms in K^^^ : the shift in R 
is tiny and easily overwhelmed by the uncertainty in ra- 
diative corrections. It seems therefore more promising to 
use the more detailed experimental information contained 
in differential distributions for this purpose. In a pioneer- 
ing study, the KTeV collaboration has attempted such an 
extraction of structure dependent contributions for K^^^ 
decays [13] by measuring the photon energy distribution, 
and their results were analyzed in detail in Ref. [20]. Here 
we provide the theoretical basis for such an analysis for 
-f^ij^, together with the chiral prediction for the outcome 
of such an extraction. 

Among the various differential rates that may poten- 
tially be investigated, the one with respect to the pho- 
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Fig. 2. Photon energy distributions from inner brcnisstrahlung 
as well as the various structure dependent terms. The notation 
dFx/dE* for the various X refers to (5.1). The normalization 

factors are Nv^^a, = 100 A^m = IQp Mk / r{Ke3)- We only cut 
on the electron-photon angle, ^e"' = 20°. 



ton energy E* is predestinated for a separation of brems- 

strahlung and structure dependent terms, as the different 
behavior of both at small E* is the defining property to 
distinguish between the two contributions. 

In the decomposition of the photon energy spectrum, 
we neglect terms that are proportional to structure depen- 
dent terms squared, and wc use the working assumption 
discussed in Sect. 3 that the structure functions can rea- 
sonably well be approximated by constants. The photon 
energy distribution thus reads 



dr dEi 



IB 



dE* 



dE* 



4 

E 

i=l 



{Vi} 



dE* 



dE* 



(5.1) 



+ o{\T^^\^,AVi,AA^. 



Here dFvi / dE* is defined to be the part of the spectrum 
proportional to V^, etc., and AVi, AAi stand for the errors 
made by the approximation of constant structure func- 
tions. 

We remind the reader that V3, A3 arc suppressed 
by TOg/M^ and are therefore essentially unobservable. 
Furthermore, we find that the distributions dFy^/dE*, 
dFA^/dE* are considerably smaller than the other struc- 
ture dependent parts of the spectrum; adding to that the 
observation that (Va), {A4) arc suppressed by two orders 
in the chiral expansion with respect to (Vi,2), (^1,2), we 
neglect these structures and only discuss the effects of Vi,2, 
A\^2 for simplicity reasons. 

The remaining parts of the photon spectrum dFx /dE* 
are shown in Fig. 2. Note that in order to put all spec- 
tra in one plot, the bremsstrahlung spectrum has been 
scaled down by two orders of magnitude. While dFiQ / dE* 
displays the expected 1/E* singularity for small photon 
energies, the structure dependent parts dFx/dE* with 



X = Vi, V2, Ai perturb the pure bremsstrahlung spec- 
trum all essentially in the same way: they rise linearly 
with E* for small energies, are bent down by phase space 
at maximum photon energies, and display a maximum be- 
tween 80 MeV and 100 MeV. Even though the strength 
of the perturbation varies, the shape is nearly the same 
for all three. Only the distribution dFA2/dE* has an ad- 
ditional node in between and is therefore suppressed. 

Similarly to the neutral kaon decay channel [20], we 
find that the bremsstrahlung spectrum is perturbed by a 
function f{E*), where 



dF, 



Vi 



dE* 



3.3 X 



dF: 



V2 



dE* 



1.3 X 



dFA, 
dE* 



(5.2) 



The information on the SD terms is contained in the ef- 
fective strength {X) that multiplies f{E*), 



dF 
dE 



dF, 



IB 



dE* 



+ {X} f{El) , 



7 7 
{X) = (Fi) +0.3(^2)+ 0.7 (Ai) . 



(5.3) 



Taking into account the chiral predictions for the various 
averaged structure functions, see Table 2, we find that 
(X) is dominated by (Vi) and {Ai), both of which con- 
tribute with comparable strength. Corrections from (V2), 
{A2) are suppressed effectively by roughly a factor of 20. 
Numerically, we predict the effective strength {X) to be 



-2.2 0(p4) 
-2.2 ±0.7 0{p^) 



(5.4) 



We therefore conclude that a significant part of the struc- 
ture dependent photon spectrum is due to the chiral 
anomaly (^i). This is different from the decay K^^^ [20], 
where Ai vanishes at O(p^), and the contribution of A2 
to the photon spectrum is comparably suppressed as in 
Fig. 2. The importance of the chiral anomaly for the ex- 
istence of various contributions to T-odd correlations in 
K't^^ was emphasized in Ref. [23], therefore it plays a sig- 
nificant role in the phenomenology of this process. 



5.2 Other distributions 

It is particularly interesting to trace back effects of the 
chiral anomaly even more clearly, i.e. we would like to 
separate the contributions of Vi and Ai that can only be 
measured as a (weighted) sum in the photon energy distri- 
bution. This turns out to be rather difficult: in many dif- 
ferential distributions, both are practically indistinguish- 
able from each other and/or from the dominant brems- 
strahlung spectrum. Two notable exceptions are shown 
in Fig. 3: the distributions with respect to the (cosines 
of) the angles between electron and pion, as well as be- 
tween neutrino and photon, show strongly different behav- 
ior in backward directions for the parts of the spectra pro- 
portional to (Vi) and (Ai). In particular, dFA^^/d cos 9*^, 
dFAi /d cos go through zero here. 
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Fig. 3. Distributions with respect to cos 9*.^, cos 5*^ from in- 
ner bremsstrahlung as well as the various structure dopcudeut 
terms. The notation dFx /d cos 61^,, dFx / d cos 8^^, is chosen 
in analogy to (5.1). The normalization factors are Nvi,Ai = 
IOOMb = W^/r(Ke3). Both spectra are shown for the stan- 
30 MeV, 61^"* =20°. 



dard cuts £™* 



A strategy for an advanced study of structure depen- 
dent terms that tries to disentangle the two most impor- 
tant SD contributions (Vi) and {Ai) on purely experimen- 
tal grounds, without input from ChPT, might therefore 
proceed along the following lines: 

1. determine the effective strength {X) w (Vi) -I- 0.7{Ai) 
from an analysis of the photon energy spectrum; 

2. use this constraint for an additional fit to the angular 
spectra dF/dcosO*^ and/or dP/d cos 9*^ in order to 
determine {Vi) and {Ai) separately. 

The (numerically) sublcading structure functions (V2), 
{A2) are, due to their smallness, even more difficult to de- 
termine. In principle, one finds a similar picture for these 
as m the case of K^!^^ [20]: the distribution dryJdE* 
peaks at lower pion energies than bremsstrahlung and the 
dominant structure dependent terms, and dFA^/d cos 9*^ 
produces a slightly enhanced variation in backward di- 



Fig. 4. Distributions with respect to cos 9'^^ from inner brems- 
strahlung as well as the various structure dependent terms. The 
notation drx/dcos6*-y is chosen in analogy to (5.1). The nor- 
malization factors are iV\A,A, = 100 Mb = W^/r{Ke3). The 
photon energy cut K"' = 30 MeV was applied. 



rections. However, these effects are expected to be much 
harder to measure even than the structure dependent 
modification of the photon energy spectrum. 

We want to briefly comment on the differential distri- 
butions dr/d cos 91^, the different contributions to which 
are displayed in Fig. 4. [Note again that the brems- 
strahlung contribution is scaled down by two orders of 
magnitude with respect to the structure dependent ones.] 
In Ref. [8], a discrepancy between data and Monte Carlo 
was found for cos 6* 



67 



-1, and the authors comment 
that this "could possibly be interpreted as a direct emis- 
sion effect". It is obvious from Fig. 4 that ChPT makes 
this claim untenable: the effects of structure dependent 
terms are too small and in particular far too flat in their 
angular distributions to produce a visible peak as seen in 
Fig. 8 of Ref. [8]. The resulting total curve IB-^SD would 
still lie within the width of the grey line for the IB contri- 
bution shown in Fig. 4. 



6 Conclusions 

In this article, we have analyzed various aspects of K'^:^^ 
decays, spinning forth previous work on K^^^. Our find- 
ings can be summarized as follows: 

1. We have constructed a decomposition of the K^^^ 
decay amplitude into an inner bremsstrahlung and 
a structure dependent part, in the absence of radia- 
tive corrections, that guarantees that the SD part has 
no non-trivial analytic properties except for unitarity 
cuts. 

2. By applying this decomposition to the chiral represen- 
tation of the K^^^ amplitude of order p** [18], we derive 
the leading chiral predictions for the structure func- 
tions. The axial structure functions that are given in 
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terms of the Wcss Zumino Wittcn anomaly arc con- 
stant, while the vector structure functions, although 
given in terms of complicated loop functions, are shown 
to be free of cuts in the physical region and therefore 
smooth and very close to constant. 

3. We consider higher-order corrections of order to the 
structure functions, in particular, we present the com- 
plete order expressions for the axial terms. Although 
cuts in the physical region due to intermediate two- 
pion states appear, their cusp behavior is shown to 
be of P-wave type and therefore suppressed, such that 
the real parts of the Ai are still smooth to a reasonable 
approximation. The size of (incomplete) higher-order 
polynomial corrections to the vector structure func- 
tions is estimated to be of natural size for chiral SU(3) 
corrections. 

4. We have investigated the stability of the prediction 
of the relative branching ratio R, normalized to the 
non-radiative width, with respect to Kf,^ form factor 
parameters. Within the range of uncertainty of the 
available experimental information on these parame- 
ters, R is shown to be remarkably stable at the per- 
mille level. Structure dependent terms reduce R by 
about 1%. Our combined prediction for experimental 

20° is 



cuts of = 30 MeV, 61^"* 



R = (0.640 ± 0.008) X 10" 



(6.1) 



where the uncertainty is dominated by unknown radia- 
tive corrections. Predictions for other cut values were 
also given in Table 6. 
5. We have discussed how to extract a linear combina- 
tion of structure dependent terms from an experimen- 
tal analysis of the differential rate dF/dE*. In contrast 
to K^^^, the axial anomaly contributes significantly 
to the perturbation of the bremsstrahlung spectrum, 
such that effects may actually be visible in K^^^. As a 
second step, we have pointed out that angular distri- 
butions of the type dF / dcos9*^ and/or dF /d cos df^^ 
may allow for a disentanglement of the two dominant 
structure dependent terms Vi and Ai. 

It would be extremely interesting and rewarding to see the 
various predictions tested by the modern high-statistics 
kaon decay experiments such as NA48/2 [11] or KEK- 
E470 [12], in particular to find miambiguous signals for 
the physics behind the structure dependent terms. 
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A Inner bremsstrahlung in K^^^ decays 

In this appendix we discuss the separation of the hadronic 
tensor V^^ into an IB and a SD part. The relevant dia- 




-yem -yhad 



a) Nti" 




b) \;t» 



Fig. A.l. Diagrams for V^u^ evaluated in the framework 
of ChPT. The hatched blobs denote one-particle irreducible 
graphs. 



grams can be grouped in two classes, displayed in Fig. A.l. 
The hatched blobs denote one-particle irreducible contri- 
butions. The diagram b) generates a pole in the variable 
u= {p — q)^, at u = M^, corresponding to the intermedi- 
ate state. We isolate the contribution of this pole by 
writing 

V = V + [2pMW') + W,MW')] , (A.1) 

where W = p — p' — q. In the following, we assume that 

this is the only singular part at g = in the tensor V^,y, 

or, in other words, that V^i, is regular at g = 0. This is the 
only assumption in the derivation of the final expression 

for the IB term. We have checked that it is true at one- 
loop order in ChPT, and we see no reason why it should 
not be correct to any order, and thus true in QCD. 
We write this regular part as 

= A=[vogi^u + vip'^q^ + V2Wfj,q^ + V3p'^W^ 

+ V4 p'^,pI + V5 W^pI + ve Wf,W,] . 

The Ward identity (2.6) generates three conditions on V^j, 
that can be written as 

vo + vi p'q + V2 qW = fx , 

Vip'q + v&qW ^ , (A.3) 
Vip'q + v^qW = 2A/+ , 

A/i = /i(i) - /.(l^') , « = (A.4) 



with 



The first equation can be solved for vq. Furthermore, we 

set 

2A/+ , _ A/i 
^5 = „„7 , Vq = -j^ + Va . (A.5) 



and find 



qW " qW 

Vip'q + V5 qW = , 
Vs p'q + qW = . 



(A.6) 
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We use the fact that the Lorcntz invariant amphtudcs Vi 
are defined for any value of the kinematic variables p'q, 
qW, and that the amplitudes are assumed to be non- 
singular at p'q = 0. It then follows that 1)5,6 are propor- 
tional to p'q, 

V5,6 = -p'qvi^a , (A. 7) 

where the sign and the numbering is chosen for conve- 
nience. Finally, we obtain 

V3A = qWv3,4 . (A.8) 

Collecting the results, wc find that V^^ can be written in 
the form displayed in (2.14), with 

(^1,^2,^3,^4) = {Vl,V2,V3,V4) . (A.9) 

(3.3) contains the explicit expression of the form factors Vi 
in the limit q = 0, t = 0, illustrating that they indeed are 
non-singular at g = at next-to-leading order in ChPT, 
as mentioned above. 

B Traces 

Here, we give the explicit expression for the sum over spins 
in |Tp in the limit where the relevant traces are evaluated 
at rrie = 0. We write 

Yl \Tf = «i f+itf + «2 f+it)Sf+ + as 5fl 

spins 

4 

+ y\{hiReVi + h\ReAi)f+{t) 

i=i (B.l) 

+ (ciReVi +cf ReAi)(5/+ 
+ 0{lmVi,lmAi,V^,AlViAi) , 

with 

5U = M^qW)-' [U{t) - U{W')] , 

(B.2) 

N = 8TraG%\Vus\^Mj, . 

For the neglected terms proportional to imaginary parts 
in the structure functions Vi, Ai, see Ref. [23]. With this 
convention for iV, the right hand side in (B.l) is dimen- 
sionless. In the limit = 0, we immediately have 

bs = bl ^ C3 = 4 = . (B.3) 
We use the abbreviations 
z pp' = a , zpq = b , z ppe — c, z ppi, — d , 
zp'q = e, zp'pe = f, zp'p^ = g , zpeq = h, 

(B.4) 

zpiyq = j, zpePv = k, zpW = l, zp'W = m, 
zqW = n, zM^ = r, z = M^^ , 

and decompose all the coefficients according to = ai a^ 
etc., where the prefactors hi, hi . . . are collected in Ta- 
ble B.l. We obtain the following expressions for the coef- 



Table B.l. Prefactors that multiply the a, , hi etc. 



ai 


A/{b'h) 


Si 


4:/{bhz) 


hi 


4,/ibhz) 






62 


4/{bhz) 


6i 


A/{bhz) 


as 


4/6^ 


64 


2/{bhz'^) 




2/{hz'^) 






Cl 


A/{bz) 




A/{bz) 






C2 


A/{bz) 




4n/{bz) 






C4 


A/ibz') 


tl 


1/z' 



ficients bi and so on: 

ai = -2abgh + 2b'^f g + 2b'^ge-b^jr 
—b^kr + 4:bcfg + 2bcge — bcjr 
-2bckr + bdhr - 2 f gh + hkr , 

0,2 = 2abghn + Ab'^fgk-2b'^ghm + 2b'^gke 
^b^hkr - h^jkr -2 b^k^r -Abcfgn 
~2bcgne + bcjnr + 2bcknr — bdhnr 
-Abfghl + 2bhklr + 4fghn-2hknr, 

03 = -U'^fgk + 2b^k^r + 4:bfgln 
—2bklnr — 2 f gn^ + kn^r , 

bi = -afhj-agh'^ + bfj + bfgh-bfke 

+2b g he — bh j r + eg he + d f he , 
62 = —ahkn + bfjk + bghk + bghn 

—bhjm + bhke — bk^e + eghn 

+d f hn — f hj I — g h^l + hkle, 
bi = —4a f g hn + 2ahknr + 4b f^gn — 4b f g ke 

+2b f gne — b f j nr — 2b f knr + 2bghme 

—bghnr — 2bg k e^ — bhkre + bjkre 

+2bk'^re + 4fghle~2hklre, 

bl = afhj-agh'^-bfj + bfgh + bfke + 2bghe 

—bhjr — bhkr + cghe — chjr — dfhe + dh^r , 
62 = ah^k — ahjk — bfhk + bfjk — bghn 

+bhjm + bhkm — bhke — bk^e + chjm 

—ehke — dh^m + dhke , 
6| = 2 f gne — fjnr + 2ghme — ghnr 

—2gke^ — hkre + j kre , 
Cl = afjn + aghn — bfjm + bfke 

—bghm + bgke — cgne — dfne, 
C2 = akn"^ — 2bfjk + bfkn — 2bghk 

+b gkn + 2bk^e — bkran — cgri^ 

—dfn^ + f j In + g hln — klne , 

C4 — 2afgn^ — ak n'^r + 4bfgke — 2bfgmn 
—2b k^r e + bkmnr — 2fglne + klnre, 

cl = —afjn + aghn + bfjm — bfke 
—bghm + bgke + bhkr — bjkr 
—cgne + cjnr + dfne — dhnr. 
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+ ^-^{CLMl + CWMl). 
The loop function Hab{x) is given by 



(C.l) 



n 7-19 1 'JabyXj'T' ^. o 



X MlMl_xSab^mlMl Ml 

.A oo_9. iOg^2 



(C.2) 



where the loop function Jab{x) is defined according to 



d) 



Jab{t) = Jab{t) - Ja6(0) 



(C.3) 



Jab{q^) 



1 f dfil 



e) 

Fig. C.l. Diagrams that contribute to the anomalous amph- 
tude A^i, at order [in a), the contribution from tt, 77 inter- 
mediate states vanishes] . Charges of the mesons running in the 
loops are not indicated. The filled vertices denote a contribu- 
tion from the anomalous Lagrangian at order p^. External line 
insertions in the tree diagram of order are not displayed. 



The other functions can also be written in relatively com- 
pact forms: 



1 



24F2 



m 



1 A Ml 

'^^Kn — 8Aktt 



= —ahk + ajk — cjm + cke + dhm — dke, 



cl = 0. 



(B.5) 



-I- 



t 



C Axial form factors at order 

In this appendix, we give the explicit formulae for the 
next-to-leading order corrections to the axial form factors 
Ai, A2, and A4, as written out formally in (3.4). The 
necessary loop diagrams for this calculation arc displayed 
in Fig. C.l. We find the following combinations of loop 
functions and counterterms: 

Si{s) = H^^is) + 3Hkk{s) + C{, 8 , 



MlMl{2AK^ + t) Ml 
A-kHAk^ °^ M2 
_ {t-2,EK^){t- Ak^) 
12-kH 



(C.4) 



M^Ml S, 



A. 



riK 



-I- 



S2{S) 



IStt^ 



8 



C2sS , 



-( 



8A'^r,K\A 
'^^VK + S^rjif 



1 , , 1 , , 167r2 
167r2 



M^MI{2Ak. - t) Ml 



K 



T2{t) = TK^{t) - Tr,K{t) + ^ C^2t t , 



{t-3SnK){t- Ak^) ' 
UttH 



Ui{u) 

U2{U) 

Xi = 2/z, - /.^ - /z^ + ^ [Cl^Ml + CIkMI) 
X2 



— « log- 



JriK{t) 

(C.5) 
(C.6) 



167r2 



19 



19 



where we have used X'ab_= M^+M^, Aab = M^-M^, and 
the two point function Jab{x) as defined in (C.3). The com- 
binations of low-energy constants occurring in (3.4) and 
(C.l) are given in Table C.l according to the numbering 
in Ref. [38]. 
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Table C.l. The coefficients from (C.l) in terms of tlie 
renormalized low-energy constants Cj^"". For example, C{s = 

— C'13'' — 5C14'' + Constants without superscript r are scale 

independent. 







Clt 


CL 


ClTi 




C2S 


CJt 


C2u 




C2K 




(jWr 











24 


-24 




















(jWr 











16 


4 











64 








(jWr 











-4 


8 











-16 


































64 


-16 































-48 


































96 


-48 





^13 


-1 


-4 


-1 


-2 


4 


-4 


-2 





-22 


2 


-16 


Oi4 


-5 


7 


-2 


-6 


-3 


-20 


8 





-4 


-8 


48 




4 


-2 


4 


-6 


-6 


16 


8 





-40 


-8 


-32 


(-^19 


1 


1 


1 


-1 


-1 


4 


2 





-2 


-2 







-1 


-1 


-4 


4 


1 


-4 


-8 





20 


8 


16 


/^iWr 


-4 


-4 


-4 


4 


4 


-16 


-8 





8 


8 





'^22 


5/2 


5/2 


1 


-1 


-5/2 


-6 


2 


4 


4 


-2 


8 


/^W 
•-^23 


-9/2 


-9/2 


-6 


6 


9/2 


-6 


-6 





12 


6 


8 



D Numerical parameters 

We denote the neutral pion and charged kaon masses with 
and Mk, respectively. In numerical evaluations, we 

use 



Mk = 493.68 MeV , M„ = M^o = 134.98 MeV 



M^± = 139.57 MeV 
= 92.4 MeV , 



me = 0.511 MeV 
Fk = 1.22 F^ . 



(D.l) 



The K^z form factor is parameterized by 



U{t) = /+(o) 



As explained in the main text, the precise values of /+(0) 
and A+ do not matter in the present context. For numeri- 
cal evaluations, we use the parameter-free one-loop result 
including isospin corrections [39] 

f^^''\Q) = 1.022 X (0) = 0.998 . (D.3) 

We stick to the low-energy constants chosen in Ref. [20] , 
Ll{Mp) = 6.3 X 10"^ , L\q{Mp) = -4.9 x IQ-^ , (D.4) 
which lead to a central value for A_(_ at one loop of 

A^;_ = 0.0294 . (D.5) 



We express the low-energy constants in the following, 
scale-independent form [18]: 



Lg = Ll{^i) - 



1 



512772 
1 

5127r' 



log 



■log 



(D.6) 
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